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A higher-order dispersive equation is introduced as a candidate for the governing equation of 
a field theory. A new class of solutions of the three-dimensional field equation are considered, 
which are not localized functions in the sense of the integrabiUty of the square of the profile over 
an infinite domain. For the new type of solutions, the gradient and/or the Hessian/Laplacian are 
square integrable. In the linear limiting case, analytical expression for the pseudolocalized solution 
is found and the method of variational approximation is applied to find the dynamics of the centers 
of the quasi-particles (QPs) corresponding to these solutions. A discrete Lagrangian can be derived 
due to the localization of the gradient and the Laplacian of the profile. The equations of motion 
of the QPs are derived from the discrete Lagrangian. The pseudomass ("wave mass") of a QP is 
defined as well as the potential of interaction. The most important trait of the new QPs is that at 
large distances, the force of attraction is proportional to the inverse square of the distance between 
the QPs. This can be considered analogous to the gravitational force in classical mechanics. 



PACS numbers: 05.45.Yv, ll.lO.-z, ll.lO.Ef 



I. INTRODUCTION 

The first half of the 20th Century oversaw the forma- 
tion of the current paradigm in physics, in which par- 
ticles and fields are the main distinct forms of matter. 
Then naturally arose the question of the interconnection 
between these two facets of the physical reality. In the 
early 1960s, Skyrme [J [2] came up with the idea that a 
localized solution of a given field equation can be consid- 
ered as a particle. Because of the formidable mathemat- 
ical difficulties in multiple spatial dimensions, the new 
idea was demonstrated in ID for the case of the sine- 
Gordon equation (sGE), for which an analytical solution 
was found in [2] for a profile composed by the superposi- 
tion of two localized shapes. The shapes interacted (scat- 
tered) in a very similar fashion as two particles would do 
if they colided with each other. In a seemingly unrelated 
research, Fermi and coworkers "^ discovered numerically 
that a virtually random initial condition for the differ- 
ence equations modeling an atomic lattice tends in the 
long term to organize into a chain of localized shapes. 
Observing that the Korteweg-de Vries equation (KdVE) 
is the limit of the difference equations modeling the lat- 
tice, Zabusky and Kruskal |4J performed similar numer- 
ical experiments for KdVE and confirmed the tendency 
of the initially "thermalized" modes to organize into lo- 
calized waves that interact as particles: the term soliton 
was introduced for this kind of wave. 

Nowadays, the subject of solitons enjoys considerable 
attention. The surveys [SI [B] give a good perspective 
of the wide scope of soliton research at the end of the 
1970s. Applications to dynamics of atomic latices are 
summarized in [7]. In [5] an excellent updated survey of 
the application of the soliton concept in elasticity can be 
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found. Several general surveys are also available [9HTT]. 
The list of the equations for which soliton solutions are 
sought is ever expanding (see, e.g., [H]). 

If an initial condition is composed by the superposition 
of two localized waves, the evolution of this wave system 
results in virtual recovery of the shapes of the two initial 
localized waves after their collision but in shifted rela- 
tive positions after outgoing from the site of interaction 
(see ini El HSl m] I among others) . The collision prop- 
erty is what justifies calling these localized waves "quasi- 
particles" (QPs). A QP solution of a fully integrable 
system is a soliton in the strict sense. For the case of non- 
fuUy-integrable system, the term QP is safer, but many 
researchers use the term "soliton" in a broader sense that 
includes also the non-fully-integrable cases. When the 
physical system is described by equation(s) that conserve 
the energy and momentum, then the same conservation 
properties are inherited by the QPs, as would be neces- 
sary for actual subatomic particles. 

What would be called a "two-soliton" solution in the 
late 1960s, was given in 1956 in [7 for the sGE, and the 
resemblance between the two-soliton solution and a spe- 
cially constructed superposition of one-solitons could be 
seen in the cited paper. In [T], the localized solutions were 
qualitatively related to mesons assuming that the meson 
field is governed by the sGE. What is more important is 
that Skyrme introduced the fundamental notion that the 
localized (in an appropriate sense) solutions can be con- 
sidered as particles of the field described by the particular 
equation under consideration. Using the techniques for 
finding analytically the two-soliton solution, an impor- 
tant analytical result was obtained in TS] T^ by extract- 
ing the positions of the centers of the individual solitons 
from the actual two-soliton profile. Thus the trajectories 
of Skyrme's "particle" (the QP in the modern terminol- 
ogy) were found and shown to bend during the inter- 
action. The last result completes, in a sense, Skyrme's 
proposal that the dynamics of the localized solution of 



the respective field equation can describe subatomic par- 
ticle interactions. The impressive analytical success here 
can be clearly attributed to the full integrability of sGE 
and to the ID nature of the solution. This line of research 
is currently actively pursued with the goal of creating a 
kind of "multi-body" formalism abased on the connec- 
tion to the field equations [TJ. Thus, Skyrme's idea to 
establish the relation between the dynamics of soliton 
solutions of a field equation and subatomic particles is a 
fruitful paradigm (see, e.g., [I8]). 

When the field equation is not fully integrable, the 
way to obtain information about the QP behavior of the 
solution is to make an ansatz that the two-soliton solu- 
tion is fairly well approximated by the linear superpo- 
sition of two one-soliton solutions with priory unknown 
trajectories, to derive a "discrete" Lagrangian with the 
trajectories as "generalized coordinates" , and to solve the 
Newtonian-like governing equations for the point dynam- 
ics of the centers of the individual solitons. This was 
proposed in [121 US] and applied to the (non-integrable) 
nonlinear Klein-Gordon equation (KGE). The trajectory 
function for the center of a QP was called a "collective 
coordinate" in [50] (see, also, [5T]) and collective variable 
(CV) [^. Nowadays, this approach is known as Varia- 
tional Approximation (VA) . The success of the simplest 
approach with one CV inspired introducing more sophis- 
ticated VA models, some involving two (and more) CVs. 
An important step in this direction was made in [53], 
where the second CV is the "width" of the QP. Such a 
choice allows one to recover, as in ^23ii24j, the relativistic 
dynamics of a single QP in the case of KGE (see [2S] for 
multiple CVs). Unfortunately, the generalization of the 
2-CV method to 2-QP interaction is not trivial. For this 
reason the 2-CV approach was tested in [52] for the sGE 
case. The integrals for the potential of interaction are un- 
surmountable without some simplifying assumption (see 
[26]). 

Despite this shortcoming, the VA is clearly a very 
promising way of establishing the wave-particle dualism 
for QPs, provided that some more pertinent models for 
the field equation are used. We would like to mention 
here that VA is one of the specific techniques to find 
an approximate solution of the field equations. The ap- 
proach of reducing the system with distributed parame- 
ters to a discrete system amounts, in general, to a Coarse 
Grain Description (CGD) '27" . From this point of view 
the VA is the most physically consistent way to create a 
CGD of a continuous system, because it retains the main 
conservation laws related to the latter. 

Before proceeding further, it is important to mention 
that recently, another avenue for further development of 
the VA has been opened. The full two-soliton solution 
shows a significant deformation of the wave profile when 
the two main solitons are close to each other. Such a 
deformation is clearly not contained in the mere super- 
position two one-soliton solutions. There is some utility 
to try to represent this deformation via the evolution of 
the widths of the QPS (as in the previous paragraph), but 



there is also another very practical approach: to consider 
what is left from the two-soliton solution when the one- 
soliton solutions are subtracted as a new particle which 
is briefly born at the site of interaction and promptly 
dies after the separation of the main lumps (the so-called 
"ghost particles" [28]). This idea was elaborated in [29] 
and offers a very pertinent perspective on the world of 
elementary particle where the birth, death, and transmu- 
tation of the particles are incessant processes. 

Apart from sGE and KGE, there are many different 
classes of dispersive nonlinear equations which possess 
localized solutions and for which the VA proved to be a 
very effective tool for approximate solutions. The nonlin- 
ear Schrodinger equation (NLS) has arisen in modeling 
the propagation of light along optical fibers [30] . Because 
of the nonlinearity, localized envelopes can appear on the 
carrier's optical frequency, which behave like QPs. The 
phenomenology is much richer for NLS, and the solitons 
for the envelope can be 'dark' or 'bright' [6]. The bright 
solitons are sometime called 'light bullets' [3TJ [35] and 
play a very important role in discretizing (digitizing) the 
information carried through the fiber. The VA was ap- 
plied to bright NLS solitons first in [33l|34] and since then 
NLS and vector NLS are the most important examples 
of the concept of QP. An excellent survey on this subject 
can be found in [6 . 

The concept of CGD of a field as an ensemble of QPs is 
clearly an important step in understanding the nature of 
wave-particle duality. Yet, there are essential differences 
between the currently known QPs and the real parti- 
cles. The most conspicuous difference is that the QPs in 
ID actually pass through each other because when the 
force is positive it remains positive up to the moment 
of collision, and then changes its sign. It is interesting 
to note that the ID space is so restrictive that even the 
two-soliton solution has a form suggesting that the in- 
dividual solitons passed through each other rather than 
scatter (see the in-depth discussion in [H]). It turns out 
j26' that only VA can shed light on what is happening, 
in the sense that the QPs scatter, but in this process be- 
cause of the restricted freedom in ID the scattered QP 
acquires the amplitude (pseudomass) of the other one, 
and the composite wave profile appears as if the QPs 
passed through each other. This effect was discovered in 
[35] and called "bouncing of mass" of the solitons for the 
case of two coupled KdVE. The repulsion or attraction of 
the QPs depends on the interaction potential defined by 
the asymptotic behavior of their tails (see [S^, among 
others). The second difference between the QPs in ID 
and the subatomic particles is that the attractive force 
decays exponentially with the distance between the QPs, 
while in reality, the attraction is given by Newton's law 
of gravitation, for which the force decays with the in- 
verse square of the distance between the centers of the 
particles. This makes it important to investigate QPs in 
dimensions higher than one, in hope of finding a poten- 
tial of interaction that is more realistic. For this purpose 
any of the known generalized wave equations containing 



dispersion should be a legitimate basic model. 

The first 3D QPs were reported for the NLS in [5X1152] . 
where the profiles of the light bullets are found numeri- 
cally for the stationary QPs and then the VA is used to 
investigate their interaction. It is well known that the 
NLS is very fertile ground for experimenting with local- 
ized shapes, because of the envelope nature of the latter. 
It is not so straightforward a matter in the traditional 
field theories, such as sGE and KGE, because a solu- 
tion of the type of the classical kink cannot be found in 
more than one dimension. Hence the simple notion of 
a lump of deformation of the field being a QP cannot 
be extended into 3D. Recently, oscillating localized 3D 
structures were found for the sGE in [371 ISH] and called 
compactons. Forcing the solution of KGE to oscillate in 
time leads to the same equation for the 3D shape as in 
the above mentioned NLS works. 

Despite their 3D nature, the localized solutions men- 
tioned in the previous paragraph exhibit the same 
exponentially decaying tails (actually, slightly super- 
exponential), and the respective force acting between 
them will decay exponentially with the increase of the 
distance between the QPs. 

The purpose of the present paper is to explore a new 
field equation and a new concept of localization of the 
solution with the aim of bringing some of the properties 
of QPs closer to the observed behavior of the objects in 
particle physics. 



II. THE THIN SHELL AS A MODEL OF A 
FIELD 

As already stated, in this paper we set to explore the 
possibility of using a novel nonlinear equation to describe 
a field. In doing this we would like to stay as close as 
possible to the spirit of quantum mechanics, i.e., as close 
as possible to the Schrodinger equation of wave mechan- 
ics, which has proven its relevance to the subject. As 
observed by Schrodinger himself [39] and elaborated in 
different papers of the present author (see [30] and the 
literature cited therein), if written for the real or imag- 
inary part of the wave function, Schrodinger's equation 
belongs to the class of Euler-BernouUi equations [41] de- 
scribing the flexural deformations of a thin plate. While, 
the Euler-BernouUi equation is the linear part of any 
shell/plate equation, the nonlinear terms depend on the 
assumptions about the nature of the deformations. In the 
theory of prolate deformations of shallow shells/plates, 
the von Karman [42^ equations are the current standard 
[43j . In the other extreme case of very steep deforma- 
tions, the shell equation is akin to the Boussinesq equa- 
tion, but the nonlinearity is different (see [30] and the 
literature cited therein). Since, the elementary particles 
have very small dimensions, the QPs that are supposed 
to describe them may exhibit even for small elevations 
deformations of order of unity and very large curvature 
(exactly the opposite to what is considered as a prolate 



deformation). We choose the model from [3D], which is 
essentially a non-prolate model, as our field equation, 
namely 



utt = Au- a3{Auf 



(3^A'^u, 
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where /3^ is the dispersion coefficient due to the stiff- 
ness of the shell (always positive) taken as the square of 
some parameter /3 for convenience. The coefficient as 
gives the relative importance of the nonlinear term. The 
membrane tension (the coefficient of the Laplacian term) 
is taken to be equal to unity, because for the case of lo- 
calized solutions over an infinite interval (no boundary 
conditions at specific points) it can be scaled out. 

Now, the spherically-symmetric, stationary solution 
satisfies 
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where —w is the curvature of the flexural deformation. 
Note that Eq. (2a) is the same as the one solved numer- 
Because of the absence of bihar- 
the result for w is the final 



ically in [3^ and [38j 
monic terms in these works, 
result. For the time being, no analytical solution is avail- 
able for this equation, and under the boundary condition 
w{r) — > for r — >■ (X) it has a solution which decays at 
infinity as e~^ /r. This solution is presented in the lower 
portion of Fig. [TJ It corresponds to the sech solution in 
ID and can be termed the "spherical sech" . 




FIG. 1: The pseudo-localized solution (upper part) and the 
localized "spherical sech" (lower part) for /3 = 1. Solid lines: 
numerically computed profile; dashed line: best fit solution 
to the linearized problem given by 2.47546(1 — e~'")/r. 

After having the solution for w, the function u, which 
can be called "the wave function", can be found by 
quadratures, namely 
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Here, one has to take C2 = in order that u decays 
at infinity. Respectively, Ci — Q ensures that there is 



no singularity in the origin of the coordinate system. As 
suggested in j40j , one can find a best-fit approximation to 
the numerical solution for w and then integrate it twice 
according to Eq. ([3]) . The result is presented in the upper 
portion of Fig. [l] It is seen that the decay of the profile of 
u at infinity is as 1/r, which is significantly slower than 
exponential decay. 

The nonlinear terms are important for shaping the pro- 
file near the origin of the coordinate system (as can be see 
in Fig. 1), but at larger distances the equation reduces 
to the following linear equation 
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which is the celebrated Euler-BernouUi equation. As 
shown in [40| this equation is akin to the Schrodinger 
equation, so interrogating Eq. Q is a good place to start. 
The Lagrangian for Eq. Q reads 
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where in the case under consideration the volume of inte- 
gration V is the entire space. Without fear of confusion 
we will use in the integral the notation — oo and oo as 
the limits of the tipple integral in what follows. 
The linearized Eq. Q has the solution (see [iP] ) 
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which is compared to the solution of the nonlinear equa- 
tion in Fig. [T] One can see that the full nonlinear solu- 
tion does not differ quantitatively from the linear solu- 
tion given by Eq. (l6]), save the fact that the latter is not 
smooth at the origin as function of the variables x,y, z. 
It decays at infinity, but is not localized in the strict 
sense, because the integral of its square diverges linearly 
at infinity with the radius of the domain of integration. 
However, the integral of the square of its gradient con- 
verges: 

{V<P)Mxdydz = An H \^^]' p'dp 
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As it should be expected, if the gradient is in L^(IR^), the 
Laplacian belongs to the same space, in the sense that 
the integral of the square of the Laplacian converges too: 
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Because of the quantitative closeness of the solution of 
the nonlinear equation to the linear one, the statements 



about the convergence and/or non-convergence of the dif- 
ferent integrals is the same for the solutions of the full and 
linearized equations. Following [ID] , we refer to this kind 
of solution as the pseudolocalized solution. The quanti- 
tative closeness to the linear solution allows us to inves- 
tigate the interaction properties of the pseudolocalized 
QPs by analytical means using the result in Eq. ([6]). 

It is interesting to observe here that in 2D, the solution 
of Eq. Q , with point symmetry, that does not exhibit a 
singularity at the origin is 



u''^\r)^Ko{r//3) + \n{r//3), 



(9) 



where Kq is the modified Bessel function of the second 
kind. Clearly, this solution does not have a pole at the 
origin, but is not a localized function, because it increases 
logarithmically to infinity at r — ^ oo. It is not even a 
pseudolocalized one in the sense of the definition of the 
present work, because the following integral diverges 

{——)\dr= / [{K'{7^))\ + 2K'{fj) + -]dr 
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This means that in 2D, pseudolocalized solutions of the 
Euler-BernouUi equation as defined here do not exist. 



III. LOCALIZED VERSUS PSEUDOLOCALIZED 
SOLUTIONS 

Eq. ([T]) belongs to the generic class of Boussinesq equa- 
tions (BEs). In different models of flows of ideal fluids 
with free surfaces or interfaces, different BEs arise. A 
generic version of Boussinesq's equation can be written 
as 
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Eq. (Ill is a generalized wave equation containing disper- 



sion in the form of a biharmonic operator of the sought 
function. The dispersive effects of the biharmonic oper- 
ator can be countered by the presence of a nonlinearity, 
and as a result a permanent wave of localized type may 
exist and propagate without change. This balance be- 
tween the nonlinearity and dispersion was first found by 
Boussinesq [44 and we shall refer to it as the "Boussinesq 
paradigm" . In different physical situations leading to the 
same generic class of equations, the nonlinear terms can 
be different. In fluid mechanics, the most general form 
of the weakly-nonlinear approximation for flows in thin 
fluid layers is discussed in [151 HS] . 

Unfortunately, for the time being we are unable to flnd 
a two-soliton solution in the 3D case of equations of the 
type of Eq. (11), just as in 2D multi-soliton solutions 
are not known, except when the equation is simplifled 
in one of the spatial directions (such as the so-called 
Kadomtsev-Petviashvili equation (KPE), see the origi- 
nal work [37] and general surveys [TUl [13] ) ■ 



The findings about the 2D locahzed solution of the 
Boussinesq equation can be summarized as follows. The 
profile of the soliton has a super exponential decay at 
infinity, namely, e~''/-y/r. This property is not robust, 
and even for very small propagation speeds, the decay at 
infinity changes to algebraic, one proportional to l/r^. 
The most disappointing property is, however, the fact 
that the steadily propagating 2D shapes are not robust 
(structurally stable) and after some time they either dis- 
sipate or blow up. This important result was first found 
by means of a special new scheme in |48j and confirmed 
quantitatively with a different numerical approximation 
P5] . All these findings rule out the strictly localized so- 
lutions as QPs. The non-robustness with respect to the 
parameter representing the phase speed of the structure 
tells us that in order to make the standing localized wave 
move, the profile behavior up to infinity has to be instan- 
taneously changed, which is hard to imagine as a physi- 
cally realistic process. 

In order to find out if pseudolocalized solutions of 
Eq. (Ill can exist we take a^ = and consider the spher- 



ical symmetric stationary solution u = G{p): 
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where Ap is the Laplace operator in spherical coordi- 



nates, already defined in Eq. 
be integrated once to obtain 



fe|. The last equation can 



G - asG^ - /?' 



J__d_ 
p^ dp 



,dG. 
dp^ 



Gi 
P 



Go 



(13) 



In order that the solution decays at infinity we set G2 — 0. 
Now, setting Ci = gives the strictly localized solution 
which has been discussed in previous paragraphs as not 
suitable for our purposes. Respectively, Gi 7^ leads to 
a pseudo-localized solution. It is important to note that 
the linearized version of the last equation has the same 
solution given in Eq. ^ for the shell-equation Q . Let us 
call the last solution "peakon" borrowing the terminology 
introduced in 50J in a different physical situation. 

A special numerical investigation of the role of the non- 
linearity for defining the shape of the pseudolocalized 
solution was recently conducted in [51]. The result is 
that the inclusion of nonlinearity does not qualitatively 
change the profile given in Eq. ([6]). The nonlinearity can 
make the profile slightly steeper at the origin, but the 
the asymptotic behavior that is essential for the inter- 
action of two QPs remains unchanged. All this means 
that, even in the case of a generic Boussinesq equation, 
the pseudolocalized solutions exist and are valid candi- 
dates for QPs. Actually, in this case, the solution $(p) is 
even a better quantitative approximation for the profile 
of the QP, hence it can be used as the basis for a coarse 
grain description of the field. 

The amplitude $(p) can be interpreted as the elevation 
of the middle 3D surface of a shell (say, along a fourth 
dimension). To an observer confined to the middle 3D 
surface it would appear as some kind of density. For this 



reason, we use a density plot in Fig. [2] In order to see 
the internal structure more clearly, we present only half 
of the ball corresponding to the main part of the peakon. 
Naturally, the latter is not exactly a ball, because it ex- 
tends to infinity but the density is smaller far from the 
center, so if we limit the number of contours to a certain 
value (sphere), the outermost contour shown is a ball. 




FIG. 2: Resected ball representing a single radial peakon for 

/3 = 1. 



IV. COORDINATE SYSTEM FOR STEADILY 
PROPAGATING LOCALIZED WAVES 

The full three-dimensional description of the interac- 
tion of QPs of the pseudolocalized type requires signif- 
icant computational resources. It is well known in ID 
that when propagating, the QPs experience contraction 
or, in some cases, elongation of the overall scale of the 
support. In 2D, the change of shape is more intricate 
(see [35]), leading to relative contraction in the direc- 
tion of motion, on the background of overall elongation 
of the scales of the solution. The first decisive simplifica- 
tion here will be to assume that when propagating, the 
QPs do not actually change their shapes, retaining the 
spherical symmetry. This is a typical assumption when 
applying the VA. 

The second important simplification is to consider QPs 
which interact only along the line connecting their cen- 
ters. Since we have already assumed that the shape re- 
tains the spherical symmetry during the interaction, this 
is not really an additional restriction on the solution. 
Even if the QPs move with arbitrary trajectories in the 
3D configurational space, the force acting between them 
will depend only on the distance between the QPs, as 
measured along the line connecting their centers. For 
this reason we restrict ourselves to the case when the QPs 
move only in one of the directions, say the z-axis. Then 



we can consider the radially symmetric cylindrical coor- 
dinates r = \J x^ + 2/^ and z. Their relation to the spher- 
ical coordinates is p := \Jr'^ + z^ and Q = arctan(z/r). 
Then 



In similar fashion from Eq. ([8]) we obtain 
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r — p cost), z = psinO. 
The cylindrically symmetric Laplace operator is given by 



_ 1^ ^ d^ 
r dr dr dz^ 
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V. COARSE-GRAIN DESCRIPTION USING 
THE FIELD OF A SINGLE QP 

For a localized wave moving in the ^-direction accord- 
ing to the law of motion z — Z{t), without changing its 
shape, one has 



4r,z,i)-A<D(v/r2 + [z-Z(i)]2), 



(15) 



where A is its amplitude. At this stage we do not consider 
the changes of the effective support of the wave (pseudo- 
Lorentzian contraction or elongation) associated with the 
motion. In other words, the speed Z{t) is small in com- 
parison with unity. Then 
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Consequently, the integral of the kinetic energy over the 
space is: 
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Note that here the relation z — Z — pcosO is used in 
manipulating the integrals. In the last line of Eq. (17 1, 
Eq. ^ is also acknowledged. Respectively, 
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where (r, z — Z) is the vector with the respective compo- 
nents. Then it follows from Eq. ^\ 
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Now for the discrete Lagrangian we obtain 
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(21) 



and the Euler-Lagrange equation for the collective vari- 
able Z(i) reads 



d r d _, 

— \m—Z\ =0, 

dr dt ^ 



(22) 



which is nothing else but Newton's law of inertia. This 
suggests that we interpret the quantity m as the mass of 
the QP concentrated in the center of the localized wave. 
Thus, we call m the pseudomass or the "wave mass" of 
a QP. The "effective" mass xnA^ is then simply propor- 
tional to the square of the amplitude of the quasi-particle 



(see Eq. (|2l])). 

Since, the solutions can have both positive and neg- 
ative amplitude, one can define also a quantity called 
pseudovolume as the integral of the curvature 



:= - 



2tt / AA<i>{p)p^sm0dpd0^4ATT 

Jo Jo 



(23) 



This means that due to the convention with the minus 
sign, in the case of positive amplitude A > 0, we have 
a positive pseudovolume. It is interesting to note that 
the magnitude of /3 does not affect the numerical value 
of the pseudovolume. The latter is just a number which 
does not depend on /3. 

The introduction of the pseudovolume allows us to un- 
derstand better the notion of anti-particle. The mass 
of an anti-quasi-particle (anti-QP) is positive, while the 
volume is negative. In the response to a force, the amti- 
QP (negative volume) will behave as matter but when 
it overlaps with a QP (positive volume) they will cancel 
over each other. 



VI. COARSE-GRAIN DESCRIPTION OF THE 
FIELD OF TWO PSEUDOLOCALIZED QPS 

The main idea of VA is to use a reasonable approxi- 
mation for the shape of the one-wave solution and then, 
from the continuous Lagrangian, derive the discrete La- 
grangian and the laws of motion for the centers of the 
QPs involved. When a CGD is attempted, the shapes 
of the QPs are taken to be the shape of a single, sta- 
tionary, propagating solitary wave. For the case under 
consideration, we stipulate that 



AL'S^iPL) + An^ipn), pk := V''^ + {z - Z^f 



(24) 



where k assumes the values L for the left QP and R for 
the right QP. Here the constraint is acknowledged that 
the QPs are supposed to move only along the axes that 
connect their centers (z-axis in the adopted notations). 
The case of two equal radial peakons (A^ — Aj^ = 1) 
with /3 = 1 is depicted in Fig.[3]for two different distances 
between their centers. 
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(b) distance between centers is 4. 

FIG. 3: Two radial peakons for Al ~ Ar = 1 and two differ- 
ent distances between their centers. /3 = 1 



In order not to obscure the main idea, we will not 



consider here the dependence of shape on the phase speed 
Ci := Zi{t). This amounts to the assumption that the 
phase speeds are small. Such an assumption is strictly 
valid only at the initial stage after the QPs are allowed 
to move under the forces that arise due to the presence 
of other QPs. 



For the time derivative of the function u from Eq. ( 24 1 , 
we find (see the one-QP case from Eq. ([l6|): 



Ut 



^Ai$'(pi 



^ Zl{z-Zl) 

PL 



-A«$'(p«)M^^^. (25) 
PR 

In some places we suppress the explicit time argument 
of Zl^h for the sake of clarity. Similarly, for the spatial 
derivatives we have 

«. = A,$'(p,)(^^^+A;,<I>'(p^)^^^^, (26) 



PL 

Ur = Al$'(pl)— + Ar^\pr) — . 
PL PR 



PR 



(27) 



A. Potential due to the Laplacian term 

Now, turning to the potential of interaction we ac- 
knowledge that there are two terms in the Lagrangian, 
containing Vm and Au, respectively. For the former we 
have 



A^ + A^ 
TT / / [{urf + {u:,f]rdrdz = n- ^ ^ 



P 



+ ALARVa{w), 



Va{w) := 27r 



oo 2 



+ [z~Zl){z-Zr) 

PL PR 

X $'(pi)$'(pH)rdrdz. (28) 



Let us now introduce a new variable t^ — z ~ Zj^ and the 
notation w — Zr — Z^. Respectively, we keep the nota- 
tion p for PL — \Jt'^ + C^ ^-iid PR = yjr'^ + (C ~ '"')^ = 
\J p^ — 2pw cos 9 + w^. Then the potential due to the 
Laplacian term is 

Vaiw) = 27: r r^^^^^^^^^<^'{p)<P'{pR)p'' sinedpde. 



"'0 



PPR 



(29) 



After performing the integration and applying some standard manipulations, we obtain 



Vaiw) 



2(3^w 



2w, 



7 (-2 + 2/3^ + 2l3w) e-*"/^ - 2(-l + f- 2/3u;)e"/^Ei(--^) + (-2 + 2/3^ - 4/3w) e'"/^Ei(-^) 



+ (-2 + 2/32 + 4(3w) e-"'/^Ei(^) + (-2 + 2/3^ + 4(3w) ln(2/3)e-™/^ + {2 ~ 2/3^ - 4(3w) \n{w)e^'"/^ 



(30) 



where 7 ~ 0.5772156649 is the Euler constant and 
Ei(u') := /_ e*/idt is the exponential integral. 

The result most closely to the interaction integrals 
which can be found in the literature is given in [53j and 
elaborated in [53], which is one of the few treatises on 
the interaction of 2D and 3D solitons. The expressions 
from the cited works differ from ours because they are 
for strictly localized QPs. Fig. |4] shows the result. The 
upper panel presents linear axes, while the lower panel 
gives the same information in log-log format. 

The most important feature of Vai^w) is that it changes 
sign. This means that it is not uniformly attractive or 
repelling. It is attractive for moderate to large distances 
between the QPs and becomes repelling when the par- 
ticles draw closer. The magnitude of the potential in 
the repelling mode is several orders of magnitude larger, 
which means that the repulsion is very strong when the 
distance between the centers of the QPs falls below the 
threshold shown in Fig. [5] The negative part of the po- 
tential, being significantly stronger than the positive part 
acts to break a possible ensemble of several QPs (an "nu- 
cleos") and, in this instance resembles the definition of 
the weak force in physics. Naturally, the much weaker 
portion of the force acting at very large distances with 
potential 1/w can be called "gravitational" because of its 
one-to-one correspondence with the respective attraction 
law in physics. 



B. Potential due to biharmonic term 

In this subsection we deal with the last term in Eq. ^ . 
To this end we observe that 






-P//3 



/32p 



(31) 



Then for the assumed profile consisting of the super- 
position of two QP's of the type defined in the Eq. ^ 
we have 



2 






Vb{w) := 20'tt I / n{p)n{pR)rdrdz = -e 

-00 Jo P 



+ ALAByb{w), 

//3 



(32) 



The last expression is obtained by using the expression 
Eq. (A9). Clearly, the potential due to the bi- harmonic 



term is positive but not strong enough to change the main 
behavior Va{w). 



C. Total potential and force 



In order to complete the description of the dynamical 
properties of the system of two QPs, we need the cross- 
mass (see the next subsection) and total potential 



V{w):^Vaiw) + Vb{w), 



(33) 



as they were defined in and Eq. (29) and Eq. (32). In 



order to save space, we will not write explicitly the ex- 
pression for the sum, because both components are well 
specified above. The total potential is shown in Fig. |4] 
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FIG. 4: The full interaction potential for /? = 0.01. 



For the total force, Fig. [5] shows its behavior and a 
comparison to different best fit functions. The result is 
presented in log-log scale, because of the great disparity 
of the magnitudes. 
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(34) 
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FIG. 5; The force of interaction between two QPs. 



D. Crossmass 

For the kinetic part of the energy of the system of two 
QPs we derive 



[utfdxdydz = -mLAlzl + -mnA^Zl 

-J ^ ^ 

+ mLRALABZLZn., (35) 



where the pseudomasses are given by Eq. (17), m^ 
mji = 1^ , and along with the pseudomass of each QP - 
have also the crossmass or the induced mass: 



mLR{w) 



{z Zl){z ^fli$'(^)$'(^^)d^dydz = 4^ 



PiPfl 
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p cos a — w cos tl 
-f \/p2 — 2pw cos 9 + w^ 



'^'{pR)psin9de\<^'ip)pdp. (36) 



The result of the integration can be found in Eq. (All I, persion parameter /3 — 0.01. 



which gives for the crossmass the following expression: 



mLR{w) 



2tt 



[e-'"/'3(4/33+4/?2iy 



+ 2/3u;2 + ^3) - 4/3^] . (37) 



Apparently, the fact the the kinetic energy Eq. (35) 



contains additional terms in the case of two QPs was 
first noticed in |55j . The term crossm,ass was coined in 
[27] for the case of interacting solitons in the case sGE 
and the pseudo- Newtonian law of interaction of two QPs 
was elaborated in full detail therein as well, including a 
numerical solution for the trajectories. Here we follow 
the line of reasoning from [27] . 

Fig. [6] depicts the dependence of the crossmass on the 
distance between the particles for the above chosen dis- 



VII. LAW OF MOTION FOR TWO 
INTERACTING QPS 

Now we move to the essence of the coarse-grain de- 
scription, namely, replacing in the original Lagrangian 
of the system Eq. ^ the sum of Eqs. (28), (32 1, and 
Eq. (35 1. As a result after acknowledging Eq. (33) we 



obtain the following coarse-grain (discrete) Lagrangian 



L 



42 I /12 
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AlArV{\Zr-Zl\) 



-nvLAlzl + ]^mRAfzl -|- xxvlrAlArZlZr (38) 
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FIG. 6: Dependence of the crossmass on the distance w be- 
tween the two QPs. 



Note that all the terms on the first line of the above 
equality are constants, so only the terms on the second 
line contribute to the Euler-Lagrange derivatives with 
respect to the functions Zj^ and Zji, namely 



^L 2 d r dZL- 






(5L 



-AlArFUZr-ZlI 
AZ 
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A.A,-^[m..^]+A?,A[m.^] 



'dt 



dt 



dt 



di 



-AlArF{\Zr- Zl\), 



(39b) 



where F is defined in Eq. ( 34 1 



The dynamical system in Eqs. (39) for the motion of 



the QPs (the centers Z]^ and Zr) differs from one gov- 
erned by Newton's law of inertia doe to the presence pres- 
ence of the crossmass terms. These terms indicate that 
the inertia of a QP is modified by the presence in its 
vicinity of another accelerating particle. As Mach put 
it [56] : every interaction in the Universe depend on all 
the surrounding matter. For reason that should be now 
obvious to the reader we will refer to Eqs. ( 39 ) as pseudo- 
Newtonian law. 

As seen in Fig. [6) the crossmass decays with the in- 
verse of the cube of the distance between the QPs, which 
means that it becomes negligible even for QPs that are 
not very far from each other and the law of acceleration 
is practically a Newtonian one. At closer distances, how- 
ever, the crossmass terms can have a significantly impact 
on the law of motion. 



We solve the equations of motion ( 39 1 numerically, in 
order to understand the mechanism of interaction be- 
tween the QPs if left to the sole action of the force acting 
between two of them. 



VIII. RESULTS AND DISCUSSION 

For definiteness we chose /3 = 0.01. This value is 
not very small, and there should be no problems with 
very thin layers in which the force varies drastically with 
slopes l/w or 1/w'*, respectively. On the other hand, this 
(3 = 0.01 is small enough to give a realistic picture of pos- 
sible application to wave mechanics where /3 is extremely 
small (of order of the square root of Plank's constant) . 

We chose for the initial positions of the left and right 
QP to be at —10 and 10, respectively. For such a separa- 
tion, the perturbation due to the the other particle at the 
site of the first particle will be small enough, hence the 
initial speeds can be set equal to zero: Zl{0) — Zr{Q). 
The first experiment is for amplitudes A^ = Ar = 1. 
Fig. [7 shows the result form the numerical integration of 



Eqs. 



39 1 for the selected set of parameters. If we find 




(a) Full scale 




(b) Zoom 

FIG. 7: Trajectories of QPs for Al = Ar = l,/3 = 0.01 and 
Zr{0) = -Zl{0) = 10, Zi{0) = Zr{0) = 0. Solid line: ZL{t); 
dashed line: ZR{t). 

the maximum of the modulus of each QP they are equal 
to each other within the precision of the computations 
(see Fig. [s]). So, QPs of the same masses (same ampli- 
tudes) scatter with the same velocities after the collision. 
The speeds of the QPs increase up to time approximately 
t = 715 after which precipitously decrease under the ac- 
tion of the negative force. After that then begin increas- 
ing again in modulus but with changed signs. This is 
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FIG. 8: Zoom of the velocities QPs for the case from Fig. [7] 
Solid line: VL{t); dashed line: VR{t). 



a classical behavior for reflecting particles and is stipu- 
lated in every model of elementary particles. Here, it 
stems from the model of QP. 

The next experiment is to consider QPs with different 
amplitudes. We select Aj^ — 1 and Aj^ — 2 (with all 
other parameters equal) and show the result in Fig. [9] 
Note that since the pseudomasses depend on the squares 
of the respective amplitudes, the mass of the left-going 
particle (initially on the right for positive abscissa) is 
four times bigger than the mass of the right-going parti- 
cle (initially on the left). Now it is clear that the more 
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FIG. 9: Trajectories of QPs for ^l = 2, A^ = 1, /3 = 0.01 and 
Zfl(O) = -2-l(0) = 10, ^i(O) = Zfl(O) = 0. Solid line: ^^(f); 
dashed line: Znit'). 

massive particle (represented by Z^) is much less affected 
by the force and is moving much slower than the particle 
Zi, (the lighter one) . For this reason the two QPs do not 
collide in the center of the coordinate system but at dis- 
tance approximately 6. This means that the left particle 
travel 16 spatial units when the right (heavier) particle 



travel just 4 spatial units. This give a rough estimate 
of the velocities as 4 times faster for the lighter particle 
whose mass is four times smaller. This is in exact inverse 
proportion to the ratio of the masses. 

In order to assess quantitatively the magnitudes of ve- 
locities of the incoming and outgoing QPs we found the 
maxima of the absolute values of the velocities. The 
graph of the velocities of the QPs is presented in this 
case in Fig. [8] Our numerical solution shows that with 
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FIG. 10: Zoom of the velocities QPs for the case from Fig. [9] 
Solid line: VL{t); dashed line: VR{t). 

high accuracy, both QPs attain maximal moduli of their 
velocities at approximately t = 363 and for this moment 
of time vl(363) = 0.209248 while vr{363) = -0.0523124. 
These maximal speeds (moduli of the velocities) can be 
assumed as the speeds at the cross section of collision. 
The ratio of the two speeds of reflection with 12 digit 
rounding is 3.99997689628 which deviates from the ex- 
pected value 4 by 5.78 • 10"^. This can be considered as 
a perfect match with the law of reflection of particles. 

Finally, we would like to address the issue of interacting 
of a QP and an anti-QP. Let set for definiteness that 
Al^I and Ar = -1. From Eq. ([34]) (see also Fig. [5]) it 
is evident that AlArF{w) < for large w, i.e. the QP 
and the anti-QP repel each other. This means that for 
an initial condition with zero initial velocity the QPs will 
gradually move farther from each other. In the region of 
intermediate w, Ai^ArF^w) > 0, which means that the 
force is attractive. Yet, to reach the stage when the QPs 
begin to attract each other, sufficiently large velocities 
are to be imposed so the inertia can carry the particle 
through the border of the region of repelling force. Fig. [5] 
shows that the change of the sign of the force around w « 
1, and it is positive down to w ~ 0.02, while for smaller 
w it is negative again. Imposing large initial velocity can 
be called "smashing" a QP into an anti-QP. 

The coupling to electromagnetic phenomena has to be 
added to the model in order to assess its relevance to 
physical reality. Yet, many of the essential features of the 
interaction are captured by this conceptually straightfor- 
ward model based on the corse-grain description. An 
enlightening discussion of the concept of anti-particle 
from the point of view of soliton theory is presented in 
[57j . The repulsion between matter (QPs) and antimat- 
ter (anti-QPs) means that a spontaneous annihilation can 
be a rare event. A special effort has to be made to smash 
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a QP and an anti-QP with very high velocity so the bar- 
rier of attractive force can be reached. When far enough, 
the QPs and anti-QPs steer clear from each other. 



IX. CONCLUSION 

In the present paper a new field equation is explored 
using the coarse-grain description (CGD), i.e., by iden- 
tifying individual elements in the field called quasi- 
particles (QPs). 

The solitary wave solutions of the new equation are 
shown to be pseudolocalized in the sense that the wave 
profile itself is not in L^(K^) (square integrable func- 
tions), but its gradient and Laplacian do belong to 
L^(M'^) space. It is shown that the analytical solution 
of the linearized equation is a very good quantitative ap- 
proximation for the shape of the solitary wave, which 
offers the possibility to identify the interaction analyti- 
cally. 

The pseudolocalized solutions are used as the shapes 
of the QPs of the new field equation and a CGD is con- 
structed via the so-called variational approximation (VA) 
in which a discrete Lagrangian is obtained for the trajec- 
tory functions of the centers of the QPs. To this end, the 
interaction potential is computed and shown to decay as 
the inverse of the distance between the QPs. Respec- 
tively, the force depends on the inverse square of the dis- 
tance between the QPs. This is radically different from 
the currently known QPs for which the interaction po- 
tential (and the force) decay exponentially with the sep- 
aration of the QPs. Thus, the first QPs that interact 
"gravitationally" are presented here. 

The force of interaction found here is not a monotone 
function of the distance between the QPs. While at mod- 
erate and large distances it is attractive, at smaller dis- 
tance it becomes repulsive. This is not an intuitive result 
from the point of view that the shape of each QP itself is 
monotone. The explanation lies in the fact that the in- 
tegrals defining the force are over a three-dimensional 
domain, and while the QPs are closer, the attraction 
from the surrounding field created by the profile can over- 
whelm the attraction from the closest elevations. The re- 
pulsive part of the force shows a conspicuous resemblance 
to the so-called weak force. The difference in our inter- 



pretation is that this is not another force, but the same 
force between the QPs which is positive (gravitational) 
at large distances, and negative (repulsive) at closer dis- 
tances with much higher magnitude. 

From the discrete Lagrangian the laws of motion for 
two interacting QPs are derived. The notion of wave 
mass (pseudomass) is introduced and an induced mass 
(crossmass) is found alongside the standard mass. The 
interpretation is that the inertia of a QP is influenced by 
the presence of an accelerating QP in its vicinity, which 
fits well in Mach's notion of all bodies in the Universe 
being connected to each other. 

Similarly to the kink and anti-kink in sGE, the QPs 
of the new equation come with positive or negative am- 
plitude. We have defined a quantity called the pseu- 
dovolume which is the integral of the curvature of the 
profile. It turns out that the the pseudo-volume is pos- 
itive when the amplitude of the respective QP is posi- 
tive. The force acting between QPs is proportional to 
the product of their amplitudes, and hence a QP and 
anti-QP repulse each other if they are at moderate to 
large distances. At the same time the masses of the two 
QPs are always positive because they are proportional to 
the squares of their amplitudes. The anti-QP reacts to 
non-gravitational forces in the same manner as the QP, 
but the difference shows only in the binary interaction 
between an anti-QP and a QP. 

The coupled equations for the laws of motion are solved 
numerically, and the dynamics are interrogated for dif- 
ferent initial velocities of the QPs. The scattering (not 
passing through each other) of QPs is confirmed, and 
information on the refiection is gathered. 

In summary, the present paper lays a claim that a sin- 
gle field model can be found whose QPs interact in a 
fashion much closer to the known properties of the inter- 
actions of real particles. 
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Appendix A: Evaluation of Integrals 

Recall that the shape of a single QP in the linear limit 
is given by Eq. Q. 



Then pl = P and pR = \/ p^ — 2pw cos Q ^ w^ . Then the 
integral for the potential V can be manipulated to the 
following 



1. Potential from the Laplacian term 



To evaluate the integral in Eq. ( 36 ) Now we can intro- 



duce again the spherical coordinates centered at z = Z^ ■ 



OO /•TT 



<^'{pW{pr)- p^sm9dpd0^ / $'(p)$'(p^) ^1- + ^^ ^- 

PPr Jo Jo ^PPr ^PPr 



p- sin 6'dpd6' 



1 



1 o^ _ 2 /-a 

<^'{p)^'{pj,) H-^ — psmedpde+ I ^^'{p)<i>'{pR)pRp sin edpde (ai) 

Jo ^ ^PR Jo Jo ^ 



We observe that 



= PW^ (PR) , -7^ ^ P^ 

96 PR dO PR 



(A2) 



which gives 



Va{w) 



2w 



Qi(p)$'(p)(p2 - w^)dp + / Q2{pW{p)dp 



(A3) 



where the terms that depend on 9 are combined in the following integrals (note that p :— cos I 



smf 



h{p)= / '^'{pr)w d0=[*(Pfl)|,.^-*(Pi?)LJ 

Jo PR 



Q2{p)^l r^'{PR)pip—de=i r 

^ Jo PR 2 Jo 



1 r 1 



e-p«//9psine'd6l = - 



$(lp 


+ w\)- 


-H\P 


~w\) 


1- 


e-P,.//3 


1 


-Pr/P 




PR 


1 


r\p+^\ 







l_e-|p+™|//3 l_f,~\p-w\/l3 



\p + w\ 



\p-w\ 



1 r I- e^P'^/'^ 
p sin ed9- - / p sin 9d0 



(A4) 



2 Jo PR 



1 /-ip+^i 



= -P-r//3- 



dp/ 



1 



2wl3 



-^^P'-^P'^^'+'"'pwdp. (A5) 



Proceeding with the manipulations of the last two integrals we obtain 



2w 



|p+^l 



[1 _ e-P-/^]dpR 



\p-A 



2w 



lp-Pfl//3ll''+^l _ A 






[\p + w\~\p-z\ + ^ (e-|P+-l//5 - e-|p--l/^)] , 



J_ f e'^Vp'-^P^P+^'p^udp = — [e-l''+"'l/'^(/? + p + z)- e-\P-'"\IH^ + Ip - u;!)] . 
2wp J_i 2w 



Thus, we can provide the final version for Eq. ( |A5[ ) 
Q2{w) = -;5^(|p + wl - |p - z|) + -^ [e-^P+-\/f^{f3 + p + z) - e-l''-"'l/^(/3 + |p - 



2w 



2w' 



w 



+ — (e-\P+^\/0 ~ e-l"-'"!/'') . (A6) 
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Now, 



1 /-"c /l „ p-|p+™l//3 \ - e-\p-w\/P-. 



'2w Jo V \p + w\ \p-w\ 

f°°r 11 1 ~i 

/ |-;^(|p+z«|-|p-2|) + — [e-l''+™l/^(^ + p + z)-e-l''-^"l/^(/3+|p~«;|)l+ — (e-l''+"'l/^-e-l''-"'l/'5)U'(p)dp 
Jq I 2w 2w '- ■' wp ^ ' J 

g-i«//j^ e™/'3Ei(-^) e-™/'3Ei(-|) e-"'/'^Ei(|) e-"'/'3Ei(-|) e"'/'^Ei(-^) 



252t 



2/3^w 2p^w 2l3^w 2/32w 

e-W/3(in(4)„in(/32)) g-™/'' ln(/?) e""'/'' ln(u;) 



2/32i(; 



4/32^ 



^2i(; 



2/32^ 



(A7) 



2. Potential from the biharmonic term 



Applying the same notations p^ = p, Pb. = \/ P^ ^ 2pit; cos 6 + w'^ we obtain 



2-k(3^ 



oo />oo 



ri(p)il(pi?)rdrdz 



27r /"^ r Z"" e-''«/^ 



/3^ Jo Jo Pfl 



sm 



Using Eq. ( A2 ) we recast the last equation to 



27r r°°l 
P^ Jo Pw 



e'P'^/^dpB. 

\p+w\ 

ztt 



27r 



P 



/32 



~p//3 



-p'dp. 



p2dp ^ f± / [g-|p-»|/0 „ g-|p+»|/^] ^-P/H^p 







/32 



li? 



Jo 



oo 



27ri(; 



(A8) 



-«)//3 



(A9) 



3. Crossmass 



For the crossmass from Eq. ( 36 ) , we evaluate the following integral 



^U \^ / ^ J 1 ^ f ^ f ^.n \ p COS^ 9 - W COS 

'^ [p)Q7n{p,w)pdp, where Qm{p,w) ^ ^ (Pfl)— ^^ ^ ^^ 

Jo V/o^ — 2pwcos6' + m;^ 



psin0d0. 



Making use of the relation Eq. (A2), we manipulate the inner integral, from the last equation as follows 



^ / ^ ;^J/ .[PCOS-U -WCOSU pSinU J^f^/ ^/ 2/, n^ 

Q,n{p.w)= / ^'{pr Y' '' ^Ae=-\<^{pn){pcos^e~wcose) 

Jo \/p2 - 2pw cos9 + w^ uj I 



^(pR)(2pcos9sme ^ wsme)de] = -(e"(''+"')/'^ - sgn(p - w)e-l''-"'l/'^ 

J w I 



(AlO) 



U'2p 



Then 



(2/32 + 2^(ii; + p) + w(«; + 2p)) el"'+''l/'^ - (2/3|u; - p\ + 2/3^ + w{w - 2p)) e-l"'-''!/'^] | 
-4/33 + e-"'/'^ (4/33 + 4^2y^ 4. 2;3w2 + w;3) 



*'(p)Qm(P,w)pdp: 



/3i(;3 



(All) 



